Confinement in the 3-dimensional Gross-Neveu model 
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We consider the TV-components 3-dimensional massive Gross-Neveu model compactified in one 
spatial direction, the system being constrained to a slab of thickness L. We derive a closed 
formula for the effective renormalized L-dependent coupling constant in the large-iV limit, using 
bag-model boundary conditions. For values of the fixed coupling constant in absence of boundaries 
A > A c ~ 19.16, we obtain ultra-violet asymptotic freedom (for L — > 0) and confinement for a 
length L (c) such that 2.07 m" 1 < L (c) < 2.82 m 1 , m being the fermionic mass. Taking for m an 
average of the masses of the quarks composing the proton, we obtain a confining legth L^f 1 which 
is comparable with an estimated proton diameter. 
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The currently accepted theory of strong interactions, 
Quantum Chromodynamics (QCD), has an intricate 
structure which makes it difficult to be directly handled. 
For this reason, simplified models have been largely em- 
ployed over the last decades in trials to obtain clues about 
the behavior of interacting quarks. Among them one of 
the most celebrated is the Gross-Neveu model [l( , which 
considers the direct four-fermions interaction. From a 
physical viewpoint, any satisfactory model describing 
strong interactions must account for the known phe- 
nomenon of confinement of the elementary constituents 
(quarks) of hadrons. Such a property has not yet been 
analytically proved, although lattice QCD gives clear in- 
dications that it indeed happens. 

In this Letter, we consider the TV-component tridimen- 
sional massive Gross-Neveu model, in the large- A limit, 
compactified along one of the spatial directions. For weak 
coupling the model presents only asymptotic freedom, as 
expected. For strong coupling, we analytically demon- 
strate the simultaneous existence of ultra-violet (short 
distance) asymptotic freedom and confinement. The con- 
fining length obtained is comparable with an estimated 
classical proton diameter. 

Although pcrturbatively non renormalizable for di- 
mensions D > 2, the massive Gross-Neveu model in 
tridimensional space has been shown to exist and has 
been constructed in Ref. Q. The idea that perturbative 
renormalizability should not be an absolute criterion to a 
model to be physically consistent comes from works done 
more than one decade ago. The results in renormalization 
theory reported in Refs. 0,0 give a theoretical founda- 
tion to support the idea that non perturbatively renor- 
malizable models can exist and have a physical mean- 
ing |5|. In fact, the Gross-Neveu model for dimensions 
greater than D = 2 has been investigated since the mid- 
dle of the seventies. Since then we know that the diffi- 



culties related to the perturbative non-renormalizability 
of the model can be surmounted, at least for D = 3, by 
performing partial summations in the bare perturbation 
series,while maintaining its initial structure (see Refs. 
El 0] ) This point has been also studied in Ref. 
The essence of the method lies in the fact that the model 
in three dimensions may acquire a naive power counting 
behavior of a just renormalizable theory, if we perform 
a partial summation, to all orders in the coupling con- 
stant, of the chains of one- loop bubbles before trying to 
rcnormalize the perturbation series directly. 

In present work, we perform such a summation for the 
large- A model and we renormalize it employing a modi- 
fied minimal subtraction scheme, where the coupling con- 
stant counter-term is a pole at the physical value of the 
argument of generalized Epstein zeta- functions |9| . This 
technique, which employs analytic and zeta-function reg- 
ularizations, has been recently used to study the bound- 
ary dependence of the coupling constant and the sponta- 
neous symmetry breaking in the compactified (ft 4 model 
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We start from the Wick-ordered massive Gross-Neveu 
Lagrangian in Euclidian metrics (in relativistic units, fi = 
c=l), 

— 7V — 9 

C =: ^x)(t^d, + m)1>{x) : +- : [i>(x)iP(x)] :, (1) 

where m is the mass, u the coupling constant, i is a 
point of R c and the 7-matrices are either 2 x 2 or 
4x4. The fermionic field of spin ^ has A components, 
tjj a (x) a — 1,2,..., A, and summation over flavor and 
spin indices are understood. In the language of Ref. Q 
a theorem has been proved stating that a family of pa- 
rameters can be exhibited such that, for A sufficiently 
large, all renormalized Schwinger functions of the type 
S2 P (rn, A, A), p=l,2,... exist. Thus we examine in the 



2 



following the large- N zero external-momenta four-point 
function, which is the basic object for our definition of the 
renormalized coupling constant. The four-point function 
at leading order in -A> is given by the sum of all chains 
of single one-loop diagrams. This sum gives for the four- 
point function at zero external momenta the formal ex- 
pression 0] 



rg>(o) = 



1 + NuZ(D) ' 



(2) 



where £(£>) is the Feynman integral corresponding to the 
single one-loop subdiagram 
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Notice that S(-D) has dimension of mass D ~ 2 and the 
coupling constant u has dimension of mass D ~ A . 

We consider the system bounded between two parallel 
planes, normal to the distance L apart from one 

another. We use Cartesian coordinates r = (x, z), where 
z is a (D— l)-dimensional vector, with corresponding mo- 
menta k = (k x , q), q being a (D — l)-dimensional vector 
in momenta space. In addition, we assume that the field 
satisfies bag-model boundary conditions 0, ^| on the 
planes x — and x = L. These boundary conditions, 
constraining the x-dependence of the field to a segment 
of length L, allow us to proceed with respect to the x 
-coordinate in a manner analogous as it is done in the 
imaginary-time formalism in field theory. This extension 
of the Matsubara procedure to describe a model compact- 
ified in a space-coordinate has been already used in the 
case of bosonic fields [ToL ITl| . Accordingly, for fermionic 
fields, the Feynman rules should be modified following 
the prescription 



dk x JL >p 



2(n + |)tt 
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/W > - — UJ r , 
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Then, the L-dependent four-point function at zero exter- 
nal momenta has the formal expression, 



rW(o,L) 



1 + NuH{D,Ly 



(5) 



where S(-D, L) is the L-dependent Feynman sum-integral 
corresponding to the single one-loop subdiagram, 



Z(D,L) = j £ 



d D ~ 1 k m 2 -(k 2 +u, 2 ) 
L ^ J (2n) D - 1 (k 2 +w 2 +m 2 ) 2 ' 



(0) 



Using the dimensionless quantities qj = kj/2irm (j = 
1,2, ...,D-1) and a = (mL)- 2 , the single bubble T,(D,L) 
can be written as 

E(D,a) = £(L>,a,s)| s=2 



4tt 2 



Ui(s, a) — U~2(s, a) — aU 3 (s, a) 



s=2 

(7) 



where 



Ui(s,a) 



n(D-l)/2-*>+2r(8-(D-l)/2) 
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and 



U2{s,a) = —- Lt/i(s - l,a), 

2 s - 1 



U 3 (s,a) 
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1 d 
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Ui{s-l >a ), 



(8) 

(9) 
(10) 



These expressions were obtained by performing the in- 
tegrals over qj in Eq. 10 using well-known dimensional 
regularization formulas [Ij] . 

Sums of the type appearing in Eq. © can be casted 
in the form 



+ OC 



E 

n=— 00 
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a(n + ^) 2 +c 2 



= 4»Zt c {v,a)-Zl (u,a), 
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(11) 



(12) 



is one of the generalized Epstein zeta-functions 0] de- 
fined for Re^ > 1/2, which possesses the following ana- 
lytic extension 



Z\ (y,a) 



faT{v) 



£(^-1/2) 
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n—1 x y \ v / 



(13) 



valid for the whole complex z^-plane. Therefore, the func- 
tions Ui(s,a) (I = 1,2,3) can be also extended to the 
whole complex s-plane, and we obtain 



Ui(s,a) 
U 2 (s,a) 

U 3 (s,a) 



h(s,D) 



T(s-^)+4W(s,a) 



(D — l)h(s- 1,D) 
2(s- l)Va 

r(s - 1- —)+4W(s - l,a) 
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where 



h{s,D) 



2 D - 2 T{s) 



(17) 



and 



W(s,a) = 2^ 

n=l 

oo 

E 



2n 



(18) 



Replacing these expressions in Eq.JTJ), we see that 
£(£>, L, s) can be written in the form 



£(£>, L, s) = ff(i).s) + G(D, L, s) 



(19) 



where H(D, s) is a term independent of a, coming from 
the first parcels between the brackets in Eas. i|14fl . I|15|l and 
(|16f) . while G(D,L, s) is the term arising from the sum 
over the parcels containing the VF-functions in these 
equations. We are using a modified minimal subtrac- 
tion scheme, where the coupling constant counter-terms 
are poles appearing at the physical value s — 2. Then 
the L -dependent correction to the coupling constant is 
proportional to the regular part of the analytical exten- 
sion of the Epstein zeta -function in the neighborhood 
of the pole at s = 2. We see from Eqs. (fT3|) and l(TB|l 
that, for s — 2 and even dimensions D > 2, H(D) is 
divergent due to the pole of the T-function. Accordingly 
this term should be subtracted to give the renormalized 
single bubble function £#(£>, L). For sake of uniformity, 
the term H(D) is also subtracted in the case of other di- 
mensions D, where no poles of T-functions are present; in 
these cases, we perform a finite renormalization. Then, 
we simply get, for any dimension D, 



n-2( s -i) 



2Dh(s- 1,D) 

s- 1 
4h(s - l,D)a d 



s - 1 
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In dimension D = 3, using the explicit form for the 
Bessel functions Kj nX j 2 {z) = y/ne~ z / '\[2z and perform- 
ing the resulting sums, we obtain (remembering that 
a = (mL)- 2 ) 



S fl (3,L) = 
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Substituting this expression into Eq. (JSJ, we obtain the 
large- N effective renormalized coupling constant, taking 
Nu = A fixed as usual, 



S (3,£) = JVrW(0,l 



A 



l + AS fl (3,L)' 



(22) 



Fig. 1 shows the renormalized single bubble (in units 
of m), S — Efl(3,L)/m, as a function of L (in units 
of m _1 ). We see from this plot and from Eq. 121|) 
that ZtmL^oo£fl(3, L) — and that £_r(3, L) diverges 
as L — > 0. This means that (see Eq. I|22|l) the fixed 
coupling constant A corresponds to the large-iV renor- 
malized effective coupling in absence of boundaries and 
that, independently of the value of A, the system presents 
ultra-violet asymptotic freedom. Moreover, X.r(3, L) has 
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FIG. 1: Renormalized single bubble (in units of m), S = 
E_r(3, L)/m, as a function of the length L (in units of m~ ). 

The interesting point is to investigate what happens 
for finite values of L, when we start from L near zero 
(in the region of asymptotic freedom). In Fig. 2, we have 
plotted the relative coupling constant y = g(3,L)/\ as 
a function of L (in units of to -1 ), taking four different 
values of A (also in units of to -1 ). For the smallest value 
considered (A = 0.2), we only see asymptotic freedom 
with g(3, L) varying from (for L = 0) up to A (for L 
x 1 Increasing the value of A (e.g. A = 10 and A = 15 in 

Fig. 2), all the curves present a peak at L = L^ax — 2.82, 
the point at which \/g(3, L) reaches its minimum value. 
This minimum value of \/g(3,L) vanishes for A = A c ~ 
19.16 (corresponding to the full line in Fig. 2) and so 
the effective coupling constant g(3, Lmax) diverges. This 
means that, for this value of A, the system is absolutely 

(c) 

confined in the slab of thickness L m ' ax - Taking greater 
values of A > A c , the minimum of \/g(3,L) becomes 

negative, g(3, L) goes to oo at a value < Lmax and 
the system would be confined in a thinner slab. However, 
as A goes to oo, L^ tends to a lower bound equal to 



in other words, for A c < A < oo the system is 
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confined in a slab of thickness lying in the interval 
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FIG. 2: Plots of the large-iV relative effective coupling con- 
stant, y = g(3, L)/X, as a function of L (in units of m _1 ) for 
four different values of the fixed coupling constant A (in units 
of m" 1 ): A = 0.2 (dashed line), A = 10 (dotted line), A = 15 
(dotted-dashed line) and A = 19.16 (full line). The dotted 
vertical lines, passing by V™ in ~ 2.07 and L$ a x — 2.82, are 
plotted as a visual guide. 

Searching for a physical interpretation of our results, 
we shall take the fixed coupling constant A = A c implying 

(c) 

that the physical confining length is L = Lmaxj the max- 
imum conceived confining length. In order to estimate 
the value of Lo, wc need to fix a value of m. We choose 
to = 3.25 x 10 n CTO -1 (~ Q.5Mev), corresponding to an 
average between typical u- and d-quark masses (recall 
that the quark content of the proton is {uud}) [15||. This 

choice gives a confining length Lq ~ 8.67 x 10 _ 12 cm. 
This value should be compared with an estimated classi- 



cal proton diameter, 2R p ~ 6.92 x 10~ 12 cm ^(|; notice 
that this value is about 10 2 times larger than the exper- 
imentally measured proton charge radius [t7|. 

Under the assumption that the Gross-Neveu model can 
be thought about as a model sharing some of the proper- 
ties of physically meaningful theories (in the sense of ex- 
perimentally testable), we can say that we have exhibited 
in a direct way the properties of ultraviolet asymptotic 
freedom and confinement, characteristics of strong inter- 
actions. This could be contrasted with the fact that in 
QCD, although asymptotic freedom is a well established 
property on theoretical grounds, confinement has been 
confirmed only by computer simulations in lattice field 
theory. Confinement has never been, to our knowledge, 
analytically proved. 

It is worth mentioning that, in Ref. |lo|. we have an- 
alyzed the dependence of the coupling constant on L for 
the large-./V <fr 4 model, the bosonic counterpart of the 
Gross-Neveu model studied here. In that case of bosons, 
no confinement is found at all, but only asymptotic free- 
dom is encountered. Comparison of the results derived in 
Ref. |l0( with those of the present work suggests that the 
effect of confinement in hadronic matter may be mainly 
due to the fermionic nature of the quarks. 

Finally, we would like to point out that the behavior 
described in this work should not be heavily dependent 
on the confining geometry. One expects to have ultravio- 
let asymptotic freedom in any case and partial or absolute 
confinement, depending on the value of A. A thorough 
study on the subject with different geometries, and in- 
cluding temperature, is in progress and will be presented 
elsewhere. 
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